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Abstract 

We investigate the constraints on the anomalous quartic ZZ^^ couphngs through the process 
pp — )• —7- pjqZX at the LHC. Taking into consideration various forward detector acceptances 
and integrated LHC luminosities, we find 95% confidence level bounds on the anomalous coupling 
parameters. We show that the bounds on these couplings are at the order of 10^^ GeV"'^ which 
are about four orders of magnitude more restricted with respect to current experimental bounds. 
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I. INTRODUCTION 



Precision measurements of gauge boson self-interactions at the LHC will be the crucial 
test of the SUl{2) x ?7y(l) gauge structure of the standard model (SM). Any deviation of the 
couplings from the SM expectations would indicate the existence of new physics. It is very 
common to investigate the new physics via effective Lagrangian approach. The theoretical 
basis of such an approach rely on the assumption that at higher energies beyond the SM, 
there is a more fundamental theory which reduces to the SM at lower energies. Hence, SM 
is assumed to be an effective low-energy theory in which heavy fields have been integrated 
out. Such a procedure is quite general and independent of the details of the model. For this 
reason this approach is sometimes called model independent analysis. 

In this paper we have analyzed genuine quartic couplings via single Z boson pro- 

duction in a 7-proton collision at the LHC Genuine quartic couplings have different origins 
than anomalous trilinear couplings. They arise from effective operators that do not induce 
any trilinear gauge boson coupling. Hence, genuine quartic couplings are free from con- 
straints on trilinear couplings. Imposing custodial SU{2)]Veak symmetry and local U{l)em 
symmetry, C and P conserving dimension 6 effective lagrangian for ZZ^^ couplings are 
given by [imj, 

C- =C^ + C, (1) 

C-o = ^aoF^.F'^^W^«W^«" (2) 

=^a,F,^F^Pw'~'^-wf (3) 

where VT*^*-' is the SU{2)weak triplet, and F^i, is the electromagnetic field strength, and Oc 
are the dimensionless anomalous coupling constants and A is the scale of new physics. The 
vertex functions generated from the effective lagrangians (2) and (3) are given respectively 
by 3 

2'Ka 



'cos2 6'h/A2 



a^iQi^u [9al3{Pl-P2) - P2aPw] (4) 



TTQ; 

-Oc [{Pl-P2){giMagul3 + 9f,l3gau) + 9a(3{Plf^P2,y + P2i,Plu) 



2 cos2 OwA-'^ 

-Pll3{gafiP2u + gauP2f,) - P2a{9piJ,Plu + QpuPlt?)] (5) 



where pi and p2 are the momenta of photons and for a convention, we assume that all the 
momenta are incoming to the vertex. 

The current best limits on couplings are provided by the OPAL Collaboration. 

These are 



0.007 GeV-2 < ^ < 0.023 GeV-^ 
-0.029 GeV-^ < ^ < 0.029 GeV-^ 



(6) 
(7) 



at 95% C.L. In 

Studying photon-induced reactions in a hadron collider is not a very new phenomena. 
The reactions such as pp — )■ p'yyp — )■ pe~^e~p p, |6|, pp — )■ p'-fyp — ?■ p fi~^fi~p 
pp — 7- p J/tfj {ip{2S))p I?! were verified experimentally by the CDF collaboration 

at the Fermilab Tevatron. These results raise interest on the potential of LHC as a photon- 
photon and photon-proton collider. Probing new physics via photon-induced reactions at 
the LHC is common in the literature. Phenomenological studies involve: supersymmetry, 
extradimensions, unparticle physics, gauge boson self-interactions, neutrino electromagnetic 
properties etc. j8l-l22|. 



II. EQUIVALENT PHOTON APPROXIMATION AND CROSS SECTIONS 

A quasireal photon emitted from one proton beam can interact with the quarks of the 
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C. Emitted quasireal photons 
25|. Their virtuality is very 



other proton and the subprocess 7g — ?■ 7gZ can occur at the 
are described by equivalent photon approximation (EPA) 
low and it is a good approximation to assume that they are on- mass-shell. Therefore to some 
extent it is possible to study 7-proton collision at the LHC. A schematic diagram describing 
this process is given in FigUl 

Any process in a 7-proton collision can be discerned from pure deep inelastic scattering 



processes by means of two experimental signatures [26| : First signature is the forward large- 
rapidity gap. Quasireal photons have a low virtuality and scattered with small angles from 
the beam pipe. Since the transverse momentum carried by a quasireal photon is small, 
photon emitting intact protons should also be scattered with small angles and exit the 
central detector without being detected. This causes a decrease in the energy deposit in 
the corresponding forward region compared to the case in which the proton remnants are 
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detected by the calorimeters. As a result of this, one of the forward regions of the central 
detector has a significant lack of energy. This defines the forward large-rapidity gap and 
usual pp deep inelastic processes can be rejected by applying a selection cut on this quantity. 
Second experimental signature is provided by the forward detectors. Forward detectors are 
capable to detect particles with a large pseudorapidity. When a photon emitting intact 
proton is scattered with a large pseudorapidity, it exceeds the pseudorapidity coverage of 
the central detectors. The detection of this intact proton by the forward detectors provides 
a distinctive signal for the 7-proton collision. 

77 collision can also be studied in the framework of EPA. In fact, these two-photon 
processes are much more studied at the LHC Two-photon processes provide a more clean 
environment with respect to 7-proton processes due to absence of the remnants of both 
proton beams. Moreover two-photon processes are generally electroweak in nature and they 
are mostly free from backgrounds coming from strong interactions. It is more probable 
for strong interactions that take part in a 7-proton process rather than in a two-photon 
process. Hence, two-photon processes generally have less backgrounds compared to 7-proton 
processes. On the other hand, for 7-proton processes energy reach and effective luminosity 
are much higher than for two-photon processes . This feature might be important 

in probing new physics especially when the energy dependences of the anomalous cross 
sections are very high. For instance, anomalous ZZ77 couplings are described by effective 
lagrangians (2) and (3) which have an energy dimension of 6. Therefore anomalous cross 
section containing the ZZ77 vertex has a higher momentum dependence than the SM cross 
section. From a simple dimensional analysis we deduce that its momentum dependence is 
higher up to a factor of p'^ than the SM cross section where p is the momentum incoming 
(or outgoing) to the vertex. Thus, new physics contribution to the cross section rapidly 
increases when the center-of-mass energy increases and the processes which have a higher 
energy reach are expected to have a high sensitivity to new physics. Another factor which 
has to be considered in a 7-proton or a 77 collision is the survival probability. It is basically 
defined as the probability of the scattered protons not to dissociate due to the secondary soft 
interactions [27]. In the 7-proton collision with quasireal photons virtuality of photons is 
very small. In the EPA that we have considered typical photon virtuality is (Q^) ~ O.OlGeV^ 



25| . Therefore proton impact parameter is much bigger than the range of strong interactions 



and proton survival probability is expected to be large. On the other hand, the survival 
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probability for a 7-proton process is usually smaller than for a two-photon process [27|. 

ATLAS and CMS collaborations have a program of forward physics with extra detectors 
located at distances of 220m and 420m from the interaction point 
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2911 ■ These forward 



detectors have a capability to detect intact scattered protons with momentum fraction loss 
in the interval ^min < < imax which is called the acceptance of the forward detectors. The 
acceptan ce p roposed by the ATLAS Forward Physics (AFP) Collaboration is 0.0015 < ^ < 



0.15 
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29| . The acceptance of the forward detectors in CMS is similar. There are also 



other scenarios with different acceptances. CMS-TOTEM forward detector scenario spans 
0.0015 < ^ < 0.5 |8,Q. 

The process 7g — )■ 7gZ takes part as a subprocess in the main reaction pp — )■ p'-yp — )■ 
P'jqZX. We consider totaly 10 subprocesses for different type of quarks and anti-quarks: 



(i) •yu — 7- •yuZ 

(ii) 'jd — )■ '-fdZ 

(iii) 7c — )■ 7cZ 

(iv) 7s — )■ •ysZ 

(v) 76 — 7- '-fbZ 



(vi) 7u — !■ juZ 

(vii) 'yd — )■ 'jdZ 
(viii) 7c — )■ ■ycZ 

(ix) 7s — i- '-fsZ 

(x) 76 — )■ 76Z 



(8) 



Each of the subprocesses is described by seven tree- level diagrams(Fig|2]). We see from Figj2] 
that one of them contains anomalous ZZ77 vertex and others are SM contributions. The 
cross section for the main process pp — p'-fp — > pyqZX can be obtained by integrating the 
cross sections for the subprocesses over the photon and quark spectra: 



cr [pp — )• P7P — !■ pyqZX) = 



XI . 



Xl > 



dxi I dx2 
'0 



dN^ 



dN„ 



a. 



E 



dz 2z 



^max 



dxi J \ dx2 



dxi ( dN^ \ ( 21 ^ - 12 



(9) 



where X\ is the fraction which represents the ratio between the scattered equivalent photon 
and initial proton energy and X2 is the momentum fraction of the proton's momentum carried 
by the quark. is the equivalent photon spectrum (see Appendix) and is the quark 
distribution function of the proton. The summations in (jH]) are performed over subprocesses 
in dH]). The second integral in ([9]) is obtained by transforming the differentials dx\dx2 into 
dzdxi with a Jacobian determinant 2z/xi where z = ^X\X2 — \/ s/s. Minv is the total 
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mass of the final particles of the subprocess 7g — )■ 'jqZ. Nq{z^ /xi) is evaluated at 
X2 = lx\. At high energies greater than proton mass, ^ ~ xi holds. Therefore it is a good 
approximation to assume that x\raax = imax and xi mm = imin in the first integral in (j9]). 
During calculations, the virtuality of the quark is taken to be = niz^ where mz is the 
mass of the Z boson. In our calculations parton distribution functions of Martin, Stirling, 
Thorne and Watt [3l| have been used. 



In Figsl3] and H] we plot the total cross section of the process pp — )■ p7p — )■ pjqZX 
as a function of anomalous couplings and ^ for the acceptances of 0.0015 < ^ < 0.5 
and 0.0015 < < 0.15. We observe from these figures that cross sections are large for 
0.0015 < ^ < 0.5 compared with 0.0015 < ^ < 0.15 as expected. We also see from these 
figures that, deviation of the anomalous cross section from its SM value is larger for the 
coupling than j^. Therefore sensitivity limits on the coupling are expected to be 
more restricted than the limits on j^. 

In all results presented in this paper we assume that center-of-mass energy of the proton- 
proton system is y/s = 14 TeV and cross sections have evaluated numerically by a computer 
code GRACE [321 . 



III. LIMITS ON THE ANOMALOUS COUPLINGS 

During statistical analysis we use two different method. We employ a simple one- 
parameter test when the number of SM events is greater than 10. On the other hand, 
we employ a Poisson distribution when the number of SM events is less than or equal 
to 10. For AFP and CMS-TOTEM scenarios SM cross sections for the main reaction 
pp — )■ pyp — )■ p'yqZX are 0.0046 pb and 0.0047 pb respectively. Hence, the number of 
SM events exceeds 10 for integrated luminosities which are equal or greater than 100/6"^. 
Therefore for AFP and CMS-TOTEM scenarios we employ both type of the statistical anal- 
ysis depending on the luminosity. Forward detectors have a capability to detect protons 
in a continuous range of ^. Therefore one can impose some cuts and choose to work in a 
subinterval of the whole acceptance region. Imposing such cuts on forward detector accep- 
tance is useful in suppressing the SM contribution. In addition to AFP and CMS-TOTEM 
scenarios we will consider 0.1 < ^ < 0.15 and 0.1 < ^ < 0.5 subintervals of the whole AFP 
and CMS-TOTEM acceptance regions. For these acceptances the number of SM events is 
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less than 10. Therefore it is very appropriate to set bounds on the couphngs using a Poisson 
distribution. 

For the acceptances of 0.0015 < ^ < 0.5 and 0.0015 < ^ < 0.15 with a high luminosity, 
analysis is performed. The function is defined by 

x^=f "^"""/" y (10) 

where a an is the cross section containing new physics effects and 5 = is the statistical 
error. The expected number of events has been calculated considering the leptonic decay 
channel of the Z boson as the signal N = S x E x asM x -^mt x BR{Z — )■ ii), where i = e~ 
or fi~ , Lint is the integrated luminosity, E is the jet reconstruction efficiency and S is the 
survival probability factor. We have taken into account a jet reconstruction efficiency of 



E = 0.6 and survival probability factor of S* = 0.7. 



his survival probability factor was 



33| . We assume that same survival 



proposed for the single W boson photoproduction 
factor is vahd for our process. ATLAS and CMS have central detectors with a pseudorapidity 
coverage |?7| < 2.5. Therefore we place a cut of |?7| < 2.5 for final state particles. Moreover, 
we also demand that the transverse momenta of the final state photon and quark are greater 
than 15 GeV. For the acceptances of 0.1 < ^ < 0.15 and 0.1 < ^ < 0.5, we employ a Poisson 
distribution. Sensitivity limits are obtained assuming the number of observed events equal 
to the SM prediction, i.e.. Nobs = S x E x ctsm x Lint x BR{Z — )■ £i). Upper limits of 



number of events N^p at the 95% CL. can be calculated from the formula 
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34| 



Nobs 

J2 PPo^sson{Nup] k) = 0.05 (11) 

A:=0 

Depending on the number of observed events, values for upper limits Nup can be found in 



Table 33.3 in Ref. 35|. In Table [T] we present number of observed events and upper limits 
of number of events for the cases in which the Poisson distribution has been used. In Table 
U the calculated Nobs values are rounded to the nearest integer. For instance, for forward 
detector acceptance of 0.1 < ^ < 0.5, Nobs = 0.69 and 1.38 for Li^t = lOO/^^ and 200/6"^ 
respectively. Both of the Nobs values have been rounded to 1. The upper limits of number 
of events Nup can be directly converted to the limits of anomalous couplings and -p-. In 
Tables HTl and UTTl we show 95% CL. sensitivity limits on the anomalous couplings ^ and 
for various integrated luminosities and forward detector acceptances of 0.0015 < ^ < 0.5, 
0.1 < ^ < 0.5, 0.0015 < ^ < 0.15 and 0.1 < ^ < 0.15. We see from Tables [TTl and [ml that 



our limits are at the order of 10~^ GeV~^ and limits on the coupling are more restricted 
than the limits on j^. 

In this paper we considered all tree-level SM contributions for the subprocess 7g — )■ 'yqZ 
(Fig|2]). These constitute major SM contributions. Any other SM contribution coming 
to this subprocess is at the loop level and can be neglected compared to tree-level SM 
contributions. The leading order background process might be the pomeron exchange. A 
pomeron emitted from one proton beam can interact with the quarks of the other proton 
and same final state can occur. But when we analyze in detail we see that this background 
process is expected to has a minor influence on sensitivity bounds. In the deep inelastic 
scattering the virtuality of the struck quark is very high. During calculations in this paper, 
the virtuality of the struck quark is taken to be = uiz^ where mz is the mass of the 
Z boson. Therefore, when a pomeron strikes a quark it probably dissociates into partons. 
These pomeron remnants can be detected by the calorimeters and background from pomeron 
exchange can be eliminated. Furthermore, survival probability for a pomeron exchange is 
considerably smaller than that for a photon exchange. Hence, even if the background from 
pomeron exchange can not be eliminated, it can not be much bigger (probably smaller) 
than the tree-level SM contributions for the photon exchange. Finally we would like to 
stress that our bounds are not very sensitive to backgrounds. For instance, if we assume 
that background cross section is 2 times bigger than the tree-level SM contributions, our 
limits with a 200/6~^ luminosity are spoiled approximately a factor of 1.5 for 0.1 < S, < 0.5 
and a factor of 1.3 for 0.0015 < ^ < 0.5. 



IV. CONCLUSIONS 



The process pp pr/p — )■ p'yqZX at the LHC with a center-of-mass energy of 14 TeV 
probes anomalous quartic ZZ'-fy couplings with a far better sensitivity than the current 
experimental bounds. It allows to improve the current sensitivity by almost four orders of 
magnitude. The potential of LHC to probe anomalous quartic ZZjj couplings was examined 
via weak boson fusion processes qq qgyy and qq — »• qq'yZ{^ 



fusion process pp — )■ P77P — ?■ pZZp 
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21 



27 



34| . In papers 



27 



and photon-photon 



3J] authors considered 



semi-leptonic decay channel of the final Z bosons ,i.e., ZZ — )■ jj where j refers to jets. 



On the other hand in papers |12L l2l| authors considered fuUy-leptonic decay channel of the 
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final Z bosons. The bounds obtained in 
obtained in papers 



12 



271, 



2l| are considerably weaker than the bounds 



34| . This probably originates from the fact that semi-leptonic decay 



channel of the final Z bosons has a large branching ratio compared to fully-leptonic decay 
channel. In our paper we have considered leptonic decay channel of final Z. Therefore it is 
more appropriate to compare our bounds with the bounds obtained in 
obtained in 



12 



21 



2l|. The bounds 



are of the same order as our bounds. Anomalous quartic 
couplings were also studied for future International Linear Collider (ILC) and its operating 
modes of 67 and 77. The limits expected to be obtained for such a machine are comparable 
with the LHC bounds S^- 



Appendix: Equivalent Photon Spectrum 



Taking into consideration the electromagnetic form factors of the proton, equivalent pho 



ton spectrum of virtuality and energy is given by the following formula 23-l25| 



a 



dE^dQ^ 



nE^Q^ 



Q 



2 

min 



A/ 



(A.l) 



where 



Q 



2 

min 



E{E - E^ 



QV 



2 

M 



gg = 0.71GeV2 



(A.2) 
(A.3) 



In the above formula, E is the energy of the incoming proton beam and is the mass of 
the proton. Fe and Fm are functions of the electric and magnetic form factors. /Xp is the 
magnetic moment of the proton. It is taken to be /i^ = 7.78. dQ^ integration in (lA.ip can 
be easily performed analytically. After integration over (5^ f[0]) takes the form of [Sl 



dN^ 



a 



dE^ TiE-y 



1 _ El. 
E 



Q 



2 

max 



Q 



2 

min 



(A.4) 



where the function ip is defined by 



^{x) = (1 + ay) 



-Mi + -) + E 



+c(l + J 



1 



In 



^ k(l + x 

k=i ^ 

1 — b + X 



1 



X 



+ 



4x(l + x)3 
^^/.(l + x) 



(A.5) 



where 



y E{E-E,y ^ Ql 



6 = 1- — ^ ^ -3.96, c = ^ 0.028 (A.6) 

Here Q^ax Q^j^ are the upper and lower bounds of the integration. The contribution 
to the integral above Q^ax ~ 2 GeV^ is negligible. Therefore during calculations we set 
QLx = 2 GeV\ 
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FIG. 1: Schematic diagram for the process pp p^p pjqZX. 




FIG. 2: Tree-level Feynman diagrams for the subprocess jq ^yqZ {q = u, d, c, s, b, u, d, c, s, b). 

TABLE I: The number of observed events Ngf^g and corresponding values for upper limits Nup at 
95% C.L. The calculated N^bs values are rounded to the nearest integer. 
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FIG. 3: Total cross section of pp jrfp p^qZX as a function of anomalous coupling ^ for 
two different forward detector acceptances stated in the figure. The center-of-mass energy of the 
proton-proton system is taken to be ^/s = 14 TeV. 



TABLE II: 95% C.L. sensitivity bounds of the couplings and for various LHC luminosities 
and forward detector acceptances of 0.0015 < ^ < 0.5 and 0.1 < ( < 0.5. The center of mass 
energy of the proton-proton system is taken to be -s/s = 14 TeV. 

0.0015 < C < 0.5 O.K ^ < 0.5 

L{fb-^) f^(xl0-6 GeV-2) ^(xlO'^ GeV-^) ;|^(xl0-6 GeV'^) ^{xlQ-^ Ge\-^) 
30 -4.0;4.0 -6.4;6.5 -2.8;2.8 -4.1;4.1 

50 -3.4;3.4 -5.6;5.7 -2.2;2.2 -3.1;3.1 

100 -2.5;2.5 -4.1;4.2 -1.9;1.9 -2.7;2.7 

200 -2.1;2.1 -3.4;3.5 -1.2;1.2 -1.7;1.7 
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FIG. 4: Total cross section of pp jrfp p^qZX as a function of anomalous coupling ^ for 
two different forward detector acceptances stated in the figure. The center-of-mass energy of the 
proton-proton system is taken to be ^/s = 14 TeV. 



TABLE III: 95% C.L. sensitivity bounds of the couplings and ^ for various LHC luminosities 
and forward detector acceptances of 0.0015 < ^ < 0.15 and 0.1 < ^ < 0.15. The center of mass 
energy of the proton-proton system is taken to be -s/s = 14 TeV. 

0.0015 < ^ < 0.15 0.1<C<0.15 

L{fb-^) f^(xl0-6 GeV-2) f^(xl0-6 GeV-^) ;|^(xl0-6 GeV'^) ^{xlQ-^ Ge\-^) 
30 -10.9;11.0 -16.2;16.2 -9.5;9.5 -13.7;13.7 

50 -9.5;9.6 -14.2;14.2 -7.3;7.3 -10.5;10.5 

100 -6.9;7.0 -10.3;10.3 -6.4;6.4 -9.2;9.2 

200 -5.8;5.9 -8.6;8.7 -4.3;4.3 -6.1;6.1 
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